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We investigate the dynamics of entanglement given by the concurrence of a two-qubit system 
in the non-Markovian setting. A quantum master equation is derived which is solved in the eigen 
basis of the system Hamiltonian for X-type initial states. A closed formula for time evolution of 
concurrence is presented for a pure state. It is shown that under the influence of dissipation that 
non-zero entanglement is created in unentangled two qubit states which decay in the same way as 
pure entangled states. We also show that under real circumstances, the decay rate of concurrence 
is strongly modified by the non-Markovianity of the evolution. 
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I. INTRODUCTION 


The important ingredient for quantum computation 
and information processing is the presence of coherent 
superpositions. A single isolated two-level system can 
be prepared in a coherent superposition of |0) and |1) 
states, and the manipulation of such states leads to new 
possibilities for storage and processing of information [ij . 
In contrast to the ideal isolated case, the interactions of 
real quantum systems with their environment lead to the 
loss of these coherent superpositions, in other words, de¬ 
coherence. However, the more realistic case would be 
manipulation of many qubits. Coherent superposition of 
such states leads to the concept of entanglement, which 
forms a precious resource for quantum computation and 
information. The fragility of entanglement is due to the 
coupling between a quantum system and its environment; 
such a coupling leads to decoherence, the process by 
which information is degraded ii- In fact, decoher¬ 
ence is one of the main obstacles for the preparation, ob¬ 
servation, and implementation of multi-qubit entangled 
states. The intensive work on quantum information and 
computing in recent years has tremendously increased the 
interest in exploring and controlling decoherence effects 
[Mil. In this work we address the problem where each of 
the two qubits are dissipatively coupled to a local bosonic 
bath; in quantum optical sense it would mean that both 
the two-level systems are subject to spontaneous emis¬ 
sion and would imply there exist relaxation between the 
excited state to ground state. Dissipation can assist the 
generation of entanglement [l^ - [l^ that can be used for 
various quantum information processing. For example, 
F. Verstraete et al 0 have shown that dissipation can 
be used as a resource for the universal quantum computa¬ 
tion without any coherent dynamics needed to implement 
it. Contrary to other methods, entanglement generation 
by dissipation does not require the preparation of a sys¬ 
tem in a particular input state and exists, in principle, 
for an arbitrary long time. These features make dissi¬ 
pative methods inherently stable against weak random 


perturbations, with the dissipative dynamics stabilizing 
the entanglement. 

The effects on system due to environment can be classi¬ 
fied into the process with memory (Non-Markovian) and 
without memory (Markovian) effects [Tsl - I^ . In case of 
Markovian processes, the environment acts as a sink for 
the system information; the system of interest loses in¬ 
formation into the environment and this lost information 
plays no role in the dynamics of the system. However, 
due to memory effects in case of non-Markovian dynam¬ 
ics, the information lost by the system during the in¬ 
teraction with the environment will return back to the 
system at later time. This makes the non-Markovian dy¬ 
namics complicated. Understanding the nature of non- 
Markovian dynamics is naturally a very important topic 
for quantum information science, where the aim is to 
control a quantum system for use in technological ap¬ 
plications (2ll - l^ . In general, three time scales in an 
open system exist to characterize non-Markovian dynam¬ 
ics: (i) the time scale of the system (ii) the time scale of 
the bath given by the bandwidth of bath spectral den¬ 
sity (iii) the mutual time scale arising from the coupling 
between the system and the bath. It is usually believed 
that non-Markovian effects strong ly rely on the relations 
among these different time scales [251 - 1^ . 

In this paper we derive a quantum master equation for 
interacting qubits with local dissipation. The equation 
is derived utilizing the completeness of the eigen basis 
of the Hamiltonian representing the interacting qubits. 
The time evolution of the density matrix turns out to be 
the sum of the time evolution corresponding to individ¬ 
ual qubits with no cross terms. Next we solve this master 
equation for X-type states with the assumption that indi¬ 
vidual baths have the same properties. The main content 
of this paper will remain however the same as for other 
kind of states and assuming different bath correlation 
functions for each bath. The different bath correlation 
functions can give rise to different time scales in the dy¬ 
namics and is treated separately. Next we identify differ¬ 
ent regimes of dynamics (Markovian and non-Markovian) 
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and show that under non-Markovian regimes of the dy¬ 
namics, there exists finite entanglement in an initially 
unentangled state. This entanglement decay in the same 
way as the pure state entanglement and we find that the 
decay rate of entanglement is strongly modified by the 
non-Markovian behavior. 

The rest of the paper is organized as follows: In section 
II, we introduce the model Hamiltonian and derive the 

S antum master equation. In the past related works [H- 
, non-interacting qubits have been considered. These 
qubits are then coupled to a common bath. However, 
in this paper, we consider qubits interacting through 
isotropic Heisenberg interaction which is a general kind 
of interaction in condensed matter physics. In section HI, 
we solve the quantum master equation in the eigen-basis 
of the system Hamiltonian for a general class of initial 
quantum states under the assumption that the bath cor¬ 
relation function decay in the same way.In section IV we 
give the decay of entanglement of certain V-type state. 
Finally we conclude in section V with the remarks of 
wider context of our results. 


and the dissipative interaction of the system with bath 
is represented by the Hamiltonian 
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Hi = 

2—1 k 

= + (4) 

2 = 1 

where Bi = Y^kdikhk- Let 0(2) = rep¬ 

resent an operator defined in interaction picture with 
respect to system and bath {Ho = Hg + Hb), we can 
therefore write H: in the interaction picture under rotat¬ 
ing wave approximation as 

2 

Hi ( 2 ) = {b)B, ( 2 ) -k d~ {t)B\ ( 2 )]. ( 5 ) 

i=l 

The time evolution of the system operators can be evalu¬ 
ated using the eigen basis of the system Hamiltonian Hg 
as 


II. MASTER EQUATION FOR LOCAL 
DISSIPATION 

In this section we will first derive master equation for 
the reduced density matrix of the system which govern 
the dynamics of the system. We consider two qubits rep¬ 
resented by spin-4 particles or two level atoms coupled 
to each other via isotropic Heisenberg interaction. The 
qubits are subject to local dissipation through a coupling 
with a bosonic bath. The Hamiltonian of the two qubit 
system is 

Hg = J~^ 2 

= + III (it + (1) 

where af" = . J represents the energy scale of the 

system. The Hamiltonian Hg can be diagonalized exactly 
i.e. Hg\'tpi) = ei|'0i) where |'0i)’s are the eigenstates of the 
Hamiltonian Hg with the eigen energies and are given 
below (with notation |0) = | f) and |1) = | j,) ): 


ei 

= J; 

ill = 

|00) 


£2 

= 0; 

ii2 = 


|oi)] 

£3 = - 

-2J; 

■03 = 

1 

t-H 

0 

|oi)] 

£4 

= J\ 

■04 = 

III). 


write the 

total 

Hamiltonian (system-kbath) 


H = 

-Hg + HB+Hj 
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util) = 'y^Pijii’Ant ® l2\i’i) exp[i(ei - ej)r] (6) 

4 

at (i) = (V'i |Li 0 crtl'^j) exp[i(ei - e^-jr] (7) 

* j=i 

where Pij = |'0i)('0j| is the projection operator PijPjk = 
Pik and Pa = I. In interaction picture, the time 
evolution of the total density matrix (system and bath) 
prit) is given by the von Neuman-Liouville equation as 


^^ = -^[HJ{t),pr{t)]. ( 8 ) 

Here we have used h = 1. We can formally integrate 
equation ([5]) and write its solution as: 

Pt( 2 ) = pr( 0 ) - i [ ds[Hi{s), pt{s)]. ( 9 ) 

Jo 

Subsituting this solution back into the commutator of 
equation (ED, we get upto second order the following 
equation: 


dprit) 

dt 


i Hi{t),pT{0) 



Hi it ), 


Hi {s),Pt (s) 


■ ( 10 ) 


where Hb is the Hamiltonian for the bath 
2 

Hb='^Y1 (3) 

i—1 k 


The solution of the above equation depends on the ini¬ 
tial conditions of total density operator. We consider an 
initially uncorrelated situation, i.e Pt(0) = Ps(0) 0 ps, 
where pg and pB are respectively the density operators 
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for the system and bath. Tracing over the bath degrees of 
freedom and assuming that trB[Hi[t)pB\ = 0, we get the 
following time non-local master equation for the reduced 
density matrix: 


dps{t) 

dt 


= - ds trslHiit), [#/(s),;Ot(s)]]. 


( 11 ) 


dpsjt) 

dt 


-i[Hs,Ps{t)] 


+ / dT[‘t>,{T)[a^ {-T)ps{t),a+] 

i=ido 

+^UT-)[a-,p,{t)a+{-T)] 


(17) 


As the bath degree of freedom are infinite so that the 
influence of the system on the bath is small in the weak 
system-bath coupling case. As a cosequence, we write 
the total density operator pt{s) = Ps{s) 0 ps + 0{Hj) 
within the second order perturbation of system-bath cou¬ 
pling [l^, [28l - l3^ . The replacement of total density ma¬ 
trix pt{s) with an uncorrelated state Ps{s) 0 ps is called 
as Born approximation. Therefore, under Born approxi¬ 
mation we write 

= -y^s trB[Hi{t), \Hi{s),ps(s) 0 ps]]. (12) 


This represents the quantum master equation in the 
Schrddinger picture. The solution of the above master 
equation depends on the type of initial states. In the 
next section we find its solution for general Al-type ini¬ 
tial states. 


III. SOLUTION OF MASTER EQUATION 

In order to obtain the dynamics of entanglement of our 
two qubit system, we assume that the qubits are initially 
prepared in an X state [s^ : 


The above equation is in a form of delayed integro- 
differential equation and is therefore a time non-local 
master equation. Replacing Ps(s) with Ps{t) in this equa¬ 
tion [i^, we get time-local master equation: 

trB[Hi{t), [Hi{s), ps{t) 0 pb]]. (13) 

Assuming the bath in the vacuum state initially, i.e 
Pb = |0)(0|; using the form of the we arrive at the 

following equation: 


^ = cMt)] + C2[Mt)]- ( 14 ) 

This forms a non-trivial result. The master equations 
contains sums of Ci for each qubit and no cross terms 
with different £,;’s. This result is the same as that for 
the non-interacting qubits. Here we have (i = 1, 2) 


A(Ps(t)) = [ds{<^i{t-s)[ai (s)ps{t),atit)] 

Jo 

+ 4>1 (i - s)[a- (t), ps (t)CT+ (s)]} (15) 
and the bath correlation function is defined as 
= {Bi{t-s)Bj)o 

= ^ |5ifepexp[-iwfe(t-s)]. (16) 

k 


■ u{0) 

0 

0 

w{0) ■ 

0 

3 : 1 ( 0 ) 

y(o) 

0 

0 

2/*(0) 

3 : 2 ( 0 ) 

0 

. w;*(0) 

0 

0 

u(0) . 


(18) 


where we have used the standard basis 
{|00), |01), |10), |11)}. Since the normalization and 
positivity of Ps(0) i.e, tr{ps{0)) = 1 and Ps(0) > 0, the 
matrix elements u,xi,X 2 ,v are non-negative parameters 
with u + xi + X 2 + V = 1, y/uv > |w|, and y/xiX 2 > |p|. 
We can use more general forms of density matrix with 
all elements non-zero, this makes the master equation 
intractable analytically. Next we express the X state 
Ps(0) in the eigen basis of Hs as 


Ps(0) = a(0)|^/>i)(V’i|+fo(0)|iA2)(V'2|+e(0)|^/>3)(V’3| 

+ d{0)\tp4){tp4\ + c(0)|V'l)(V'4| + C*(0)|V’4)(V’l| 
+ li(0)|V'2)(V'3| + ^*(O)|V'3)('02| (19) 

where the various parameters of the density operator in 
the eigen basis of Hg are related to the parameters in the 
standard basis in the following way: a(0) = u(0), 6(0) = 
i[xi(0)-h 3 : 2 ( 0 )-hp(0)-fy*(0)], e(0) = i[a:i(0)-h 0 : 2 ( 0 )- 
2/(0) - 2/^(0)], /i(0) = i[a:i(0) - 3 : 2 ( 0 ) - 2/(0) + 2/^0)], 

d{0) = w(0), c(0) = w{0). Next we see that the form of 
the density matrix is invariant during the time evolution 
generated by the quantum master equation. Therefore 
we can the density matrix at time t as 

Ps{t) = a(t)|V’i)(V’i|+^(l)l^2)(^2| + e(t)|V'3)(^3| 

+ C^(l)l'04)('04| -hc(t)|V'l)(V'4| -Fc*(t)|V'4)('0l| 

+ h{t)\tjj2){tp3\ + h*{t)\'tp3){tj}2\- (20) 


Next we revert back to the Schrddinger picture with a In order to find the time evolution equations of the vari- 
change in variable t = t — s, we write ous parameters involved in equation (I20p we assume that 











the bath correlation functions have the same form 


1 


( 21 ) 

where A is the spectral width of the bath, is related 
to the microscopic system-bath coupling constant. It de¬ 
fines the relaxation time scale over which the state of 
the system: t/j ~ F”^. It can be shown to be related 
to the Markovian decay rate Fm in Markovian Limit of 
flat spectrum. This form of correlation function corre¬ 
sponds to the Lorenztian spectral density of the bath 
[l5|. Assuming that Fi = F 2 = Fm for simplicity, we 
substitute the Ps(t) as in equation (12011 in the quantum 
master equation (HU and obtain the time dependence of 
the parameters as 


a{t) 

1 

II 


(22) 

d{t) 

= d{0) + f dz[r]{z 

'.)b{z) + Y,{z 

)e{z)] (23) 

c{t) 



(24) 

hit) 



(25) 



= V{t)a{t) 

(26) 



= T,{t)a{t) 

(27) 

m 

= F+(t) + F_(t); 

r+(0 = 

i/gdz S(z); 


T-W = + 'S'-W; S2{t) = 

2Jt -I- Si{t) and the explicit forms of these functions are 
given in the appendix A. 


IV. DECAY OF ENTANGLEMENT 





FIG. 1: Decay of entanglement as measured by concurrence 
C{t) with time at different values of the parameter Q and R. 
Here we have used 9 = 7r/2. Plots for different Q values at 
(a) R = 100 (h) R = 1 (c) R = 0.01 and (d) R = 0.001. 
Plots (a) and (b) are in units of Markovian decay rate Pm, 
i.e T = Pm!, while plots (c) and (d) are in units of rescaled 
decay rate i.e t' = 


pure entangled state I'L) = cos||01) -|- sin||10). This 
state has initial entanglement C(0) = sind and at time t 
we write with the help of above results 


C{t) = 2 max{0, |2/(t)|} (29) 


In this section we study entan glem ent of a two qubit 
system by means of concurrence |34| . For a density ma¬ 
trix p, concurrence is defined as C = max{0, y/FT” — 
where ri, r 2 , ra and are the eigen values 
of matrix R in the descending order. The matrix R is 
defined as i? = p{<7\ ® cr|)p*(crj' (g) cr|) and p* represents 
complex conjugation of p in standard basis. For X state 
in the standard basis we write concurrence as [s^ 


C{t)=2 max{0,\w(t)\ -\/xi{t)x 2 (t), \y{t)\ -^/u{t)v{t)} 

(28) 


where we have u(t) = a{t), w(t) = c(t), Xi{t) = 
i[6(f) -I- h{t) + h*{t) + e{t)], X 2 {t) = ^[blt) - h{t) - 
h*{t) + e(t)], y{t) = i[6(t) - h{t) + h*{t) - e(t)], y*{t) = 
i[6(t) -|- h{t) — h*{t) — e(<)], v{t) = d{t). Next we use 
these results to investigate the decay of entanglement in 
some specific cases. First We consider the decay of the 


or 


c(f) = 2\y{t)\ = _ i_TMe-r.«)(3„, 


This forms an important result. It shows that even 
though we have initially an unentangled state C(0) = 0, 
we still have entanglement at later time t. This can be 
attributed to the dissipative interaction between the sys¬ 
tem and the bath. Let us suppose d = tt, which corre¬ 
sponds to 110) state; the effect of dissipative interaction 
( Hi{t)\10) = i?j(t)|ll)-|-i? 2 (i)| 00 ) ) results in an entan¬ 
gled state. 

Next we analyze the Markovian and non-Markovian 
regimes of the dynamics, for that we define the following 
parameters: r = Q = j, R = Therefore, 

using this parametrization we have 
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r+(<) 


r-W 


2[1 + 9Q2] 

3Q 

i?[l + 9Q2]2 


[1 — 9(5^] [1 — e cos{3QRt)] 

2i?[l + 9g2]2 

sm{3QRT) (31) 


2[1 + Q2] 

Q 

i?[l + Q2]2 


[1 — g^][l — e cos(gi?r)] 
2R[l + Q 2]2 


e sin(gi?T). 


(32) 


exponential decay in units of the rescaled time. Next 
we see that before reaching the limiting behavior of ex¬ 
ponential decay in rescaled time (I35L we observe some 
oscillatory behavior Figure [Ijc)-(d). The deviation from 
an exponential decay can be attributed to the memory 
effects developed in the two qubit system. This occurs 
clearly due to the second terms in r_(t) and r+(t). For 
g ^ 1 we may approximate this as 

1 _ cos(i?g3r')l 

- w - ‘ ’ 



In order to understand how the Markovian limit is ob¬ 
tained from the above expressions, we plot in Fig. [TJa)- 
(b) C(r) for 0 = 7r/2 with respect to dimensionless pa¬ 
rameter T for i? = 100 and i? = 1 at different values 
of g. We observe that the Markovian curve is recov¬ 
ered for R » 1 with g = 0. We can understand this 
behaviour of C{t) by looking at the different paramters 
involved. The typical time scale over which the system 
of two qubits changes is Ts ^ 1/J and the time scale over 
which the bath changes is tb ^ 1 /A while relaxation time 
scale for each qubit would be given by ~ I/Fm- It 
means i? >> 1 implies tb « tb and g < 1 implies 
tr < Ts- Thus physically R >> 1 and g < 1 would im¬ 
ply that the system evolves over a large time compared 
to very fast bath dynamics. Therefore Markovian regime 
corresponds to i? >> 1 and g < 1 and we have 


r+(T') « — 
^ 18 


t' + 


1 _ g-RQ^r' cos{3RQ^t') 


RQ^ 


.(37) 


In order for the oscillatory term to be visible, we re¬ 
quire the exponential decay term in (I36l) - (l37p to be not 
too fast giving RQ^ < 1 , but simultaneously the oscil¬ 
lation frequency should be faster than the overall decay 
envelope RQ^ > RQ^ > 1. The strongest oscillations 
therefore occur when RQ^ ~ 1 , which agree with the nu¬ 
merical plots in Figure [IJc)-(d). The deviation from an 
exponential decay can be attributed to the memory ef¬ 
fects developed initially, typical of non-Markovian behav¬ 
ior. The criterion for the strongest oscillatory behavior 
are satisfied when all the characteristic time scales are all 
approximately the same \.e tr ^ Ts ^ tb- 


r+(t) = r_(t) = ^ (33) 

and therefore we get standard Markovian limit: 

C(t) =C(0)e"5rMt_ (34) 

We observe that under the Markovian limit an initially 
unentangled state remains unentangled always. In situa¬ 
tions where the spectral width A of the bath is narrower 
than the energy scale J involved for the system implying 
g >> 1. This would mean tr « Ts- In Fig. [IKb), 
we observe for i? = 1 that as Q increases from 0, there is 
larger devaition from the Markovian dynamics of the con¬ 
currence C{t). The general trend is similar for all values 
of R as can be seen on comparison to Fig. [IKa). These 
observations suggest Markovian regime is in fact oppo¬ 
site to the regime i? < 1 and g >> 1, which we call as 
non-Markovian regime. The term that is responsible for 
this larger deviation can be attributed to first terms of 
F_(t) and r_|_(t) containing in the denominator. For 
g^ >> 1 suggest defining another time scale 


The decay of entanglement defined by C{t) at various 
values of Q for non-Markovian regime i? < 1 in terms of 
rescaled time r' = rjQ'^ is shown in the FiglTJc)-(d). We 
see that for large Q the concurrence C(r') coincides with 


V. CONCLUSIONS 

In conclusion, we have derived a quantum master equa¬ 
tion for system of two interacting qubits under the influ¬ 
ence of local dissipation. Using the assumption that the 
correlation functions have the same form for each of the 
baths, the solution of master equation is found for the 
general X-type state. The time dependence of concur¬ 
rence, a measure of entanglement is studied for a pure 
entangled state Idt) = cos|| 01 ) -|-sin|| 10 ) (a special 
case of X-type state) under both Markovian and non- 
Markovian regimes of the dynamics. It is found that for 
finite time evolution an unentangled state can go to an 
entangled state in contrast to the Markovian case where 
the unentangled state remains unentangled always. By 
identifying the parameter space, we have found that our 
results reduce to standard Markovian decay rate which 
in general is not a physically relevant regime . In the 
physically relevant regime with narrower spectral width 
as compared to J, the decay rate is better approximated 
by F(t) = which is the standard Markovian decay 
rate divided by the g^, which can be quite large in prac¬ 
tice. 

Next we compare our work with the several other works 
that studied non-Markovian dynamics of entanglement. 
Taking the example of Ref. [I 4 I , the authors derive the 
non-Markovian decay of the entanglement for the pure 
state I'k) = cos||01) -|-sin||10). The time dependence 
of concurrence is given by: 
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Appendix A 


C{t) = max{0, C'(0)G(i)} 

where 


(38) 


In this appendix we write the explicit forms of the 
various functions used in main text. 


G{t) = 


- ^-At/2 


cosh(^5) + i sinh(^i5) 
2 0 2 


(39) 


and (5 = y 1 — and Fm is the Markovian decay rate. 
Our result Eqn. (1501) is more general than above re¬ 
sult. The results in these works do not tell about 

the amount of entanglement and its decay that would be 
present in an entangled state generated from dissipation. 
To see more clearly the behavior, let us examine this in 
two limiting cases 

Weak coupling limit Tm A; This regime corresponds 
to a weak coupling regime or a very broad coupling to 
many frequency modes, which gives Markovian behavior. 
Here (5 « and the decay function Git) gives purely 

exponential behavior. To first order the decay function 
may be approximated as 

G{t) ^ (40) 

which is nothing but standard Markovian spontaneous 
decay. 

Strong coupling limit Tm ^ A: The reverse regime 
is when the linewidth of the bath is extremely narrow, 
which gives rise to strongly non-Markovian behavior. 

Here we may approximate (5 = and 


? 7 (t) = [A(1 — e'^* cos Jt) + Je ^‘sinJt] (Al) 

T,{t) = [A(1 — cos3Jt)-I-3Je“'^*sin3Jt] 

(A2) 


r_(t) = 


dz rj{z) 


(A3) 


JQ 

FmA 


2[A2 + J2] 


At — —^[1 —e ^‘cosJt)] 


A2-H J2 
2XJ 


A2-h J2 


e sin Jt 


r+(t) = ^J^dzX:{z) 


0 

FmA 


2[A2 -h9J2] 


At — 


A^ - 9 J2 
A2 -t 9 J2 
6AJ 


(A4) 
(A5) 

[1 — cos3Jt)] 


-Ai 


A2 -h 9 J2 


sin3Jt 


(A6) 


G(t) = 


cos 




(4ip_(t) 


w hich cor responds to damped oscillations at frequency 
a/AFm/ 2 and a decay envelope with rate A. Thus we see 
that in both the cases the previous results do not yield 
the scaling factor as derived in Eq. (35). S+{t) 

The current result would be important for applications 
where spontaneous emission is a serious drawback of us¬ 
ing excited states, such as for quantum information pro¬ 
cessors, quantum simulators, and quantum metrological 
applications. 


FmA 


2[A2 + J2] _ 


FmA 

2[A2 -h9J2] 


Jt 


2J 


3Jt 


A2-t J2 

A2- 

■A2 + J2® 

6J 


[1 — e ^‘cosJt)] 


-Ai 


sin Jt 


(A7) 


A2 -h 9 J2 
A^ - 9 


[1 — e ■^*cos3Jt)] 


-At 


A2 -h 9 J2 


sin3Jt 


(AS) 


[1] M. A. Nielsen and I. Chuang, Quantum Computation and 
Quantum Communication ( Cambridge University Press 
) 2000 . 

[2] M. Schlosshauer, Rev. Mod. Phys. 76, 1267 2005 

[3] W. H. Zurek, Rev. Mod. Phys. 75, 715(2003) . 

[4] J. T. Barreiro, P. Schindler , O. Giihne ,T. Monz , M. 
Chwalla , C. F. Roos , M. Hennrich and R. Blatt, Nature 
Phys. 6, 943 (2010). 

[5] S. Schneider and G. J. Milbunr, Phys. Rev. A 57, 3748 
(1998). 


[6] Q. A. Turchette, C. J. Myatt, B. E. King , C. A. Sack- 
ett, D. Kielpinski , W. M. Itano , C. Monroe and D. J. 
Wineland, Phys. Rev. A, 62, 053807 (2000). 

[7] C. J. Myatt, B. E. King , Q. A. Turchette, C. A. Sack- 
ett, D. Kielpinski , W. M. Itano , C. Monroe and D. J. 
Wineland Nature, 403, 269 (2000). 

[8] E. Knill , R. Lafiamme and W. H. Zurek, Science 279, 
342 (1998). 

[9] S. Diehl , A. Micheli, A. Kantian, B. Kraus, H. Buechler, 
P. Zoller, Nature Phys. 4, 878 (2008). 






































7 


[10] F.Verstraete, M. M.Wolf, and J. I.Cirac, Nat. Pliys.5, 
633 (2009). 

[11] H. Weimer , M. Muller , I. Lesanovsky , P. Zoller and 
H.P. Biichler, Nature Phys. 6, 382 (2010). 

[12] M. B. Plenio et al., Phys. Rev. A 59, 2468 (1999). 

[13] Beige et al., J. Mod. Opt. 47, 2583 (2000). 

[14] P. Horodecki, Phys. Rev. A 63, 022108 (2001). 

[15] H. P. Beuer and F. Petruccione, The Theory of Open 
Quantum systems (Oxford, New York: Oxford University 
Press) (2005). 

[16] R. Lo Franco, B. Bellomo, S. Maniscalco, and G. Com- 
pagno, Int. J. Mod. Phys. B 27, 1345053 (2013). 

[17] B. Bellomo, R. Lo Franco, and G. Gompagno, Phys. Rev. 
Lett. 99, 160502 (2007). 

[18] B. Bellomo, R. Lo Franco, and G. Gompagno, Phys. Rev. 
A 77, 032342 (2008). 

[19] K. M. Fonseca Romero and R. Lo Franco, Phys. Scr. 86, 
065004 (2012). 

[20] M. Q. Lone, and T. Byrnes, Phys. Rev. A 92, 011401(R), 
(2015). 

[21] B. M. Terhal and G.Burkard, Phys. Rev. A 71, 012336 
(2005). 

[22] M. Ban, J. Phys. A: Math. Gen. 39, 1927 (2006). 

[23] M. M. Wolf, J. Eisert, T. S. Cubitt, and J. 1. Cirac, Phys. 


Rev. Lett. 101, 150402 (2008). 

[24] F. Pastawski, L. Clemente and J. I. Cirac, Phys. Rev. A 
83, 012304 (2011). 

[25] P. Haikka and S. Maniscalco, Phys. Rev. A 81, 052103 

( 2010 ). 

[26] W. M. Zhang, P. Y. Lo, H. N. Xiong, M. W. Y. Tu, and 
F. Nori, Phys. Rev. Lett. 109, 170402 (2012). 

[27] M. A. Cirone, G. De. Chiara, G. M. Palma and A. Recati, 
New J. Phys. 11 103055 (2009). 

[28] H. J. Carmichael, Statistical Methods in Quantum Optics 
I (Berlin: Springer-Verlag) (2008). 

[29] H. P. Brener, B. Kappler and F. Petruccione, Phys. Rev. 
A 59, 1633 (1999). 

[30] H. P. Breuer, B. Kappler and F. Petruccione , Ann. Phys. 
(NY) 291, 36 (2001). 

[31] M. Schroder , U. Kleinekathofer and M. Schreiber, J. 
Chem. Phys. 124 084903 (2006). 

[32] E. Ferraro ,M. Scala, R. Migliore and A. Napoli, Phys. 
Rev. A 80, 042112 (2009). 

[33] T. Yu and J. H. Eberly, Quantum Inf. Comput. 7, 459 
(2007) 

[34] W.K.Woofers, Phys. Rev.Lett. 80, 2245 (1998). 


